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The late-time tail behavior of massive Dirac fields is investigated in the Schwarzschild black-hole 
geometry and the result is compared with that of the massive scalar fields. It is shown that in the 
intermediate late times there are three kinds of differences between the massive Dirac and scalar 
fields, (I) the asymptotic behavior of massive Dirac fields is dominated by a decaying tail without 
any oscillation, but the massive scalar field by a oscillatory inverse power-law decaying tail, (II) the 
dumping exponent for the massive Dirac field depends not only on the multiple number of the wave 
mode but also on the mass of the Dirac field, while that for the massive scalar field depends on 
the multiple number only, and (III) the decay of the massive Dirac field is slower than that of the 
massive scalar field. 

<^ , PACS numbers: 03.65.Pm, 04.30.Nk, 04.70.Bw, 97.60.Lf 
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(N ■ I. INTRODUCTION 

> 

' The dynamical physical mechanism responsible for the relaxation of perturbation fields outside a black hole and the 
decay rates of the various perturbations have been extensively studied since Wheeler introduced the no-hair 

theorem in the early 1970s The massless neutral external perturbations were first studied by Price, and it 

was found that the late-time behavior for a fixed r is dominated by the factor £~( 2 '+ 3 ) for each multiple moment I 
0. The massless charged scalar field was studied in Refs. and the conclusion was that a charged hair decay 

slower than a neutral one, i.e., the charged scalar hair outside a charged black hole is dominated by a i~( 2i + 2 ) tail. 
O The massless late-time tail for the gravitational, electromagnetic, neutrino and scalar perturbations had also been 
considered in the case of the Kerr black holes in Refs. On the other hand, many authors found that the 

analysis of massive fields is also physically important since massive fields can cause interesting phenomena which are 
qualitatively different from the massless case. The evolution of a massive scalar field in the Schwarzschild background 
was analyzed by Starobinskii and Novikov 0, and they found that, because of the mass term, there are poles in the 
complex plane closer to the real axis than in the massless case, which leads to inverse power-law behavior with smaller 
indices than the massless case. Hod and Piran [8( pointed out that, if the field mass /i is small, namely [iM <C 1, the 
oscillatory inverse power-law behavior 
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dominates as the intermediate late-time tails in the Reissner-Nordstrom background. We recently investigated the 
late-time tails of the massless and the self-interacting (massive) scalar fields in a stationary axisymmetric Einstein- 
Maxwell-dilaton-axion black-hole geometry, and found that the dumping exponents is independent of the rotation 
parameter and the dilaton of the black hole. 

Although much attention has been paid to the study of the late-time behaviors of the scalar, gravitational, elec- 
tromagnetic in static and stationary black-hole backgrounds, however, to my best knowledge, at the moment the 
late-time evolution of the Dirac fields has not been investigated. The aim of this paper is to study the intermediate 
late-time tail behavior of the massive Dirac fields in the Schwarzschild black-hole background and to see whether or 
not special properties exist in this case. 

The plan of the paper is as follows. In Sec. 2 the decoupled massive Dirac equations in the Schwarzschild spacetime 
are presented. In Sec. 3 the black-hole Green's function is introduced by using the spectral decomposition method 
|22|. In Sec. 4 the intermediate late-time evolution of the Dirac massive fields in the Schwarzschild background 
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is investigated. The section V is devoted to a summary and conclusion. In the appendix we study whether the 
conclusions might change if the tortoise coordinate is defined in the conventional way. 

II. DIRAC EQUATION IN THE SCHWARZSCHILD SPACETIME 

Dirac equation in a general background spacetime can be expressed as |23| 

[ 7 a e^(^ + r Al )+ M ]* = 0, (2.1) 

where /i is the mass of the Dirac field, r y a is the Dirac matrix, is the inverse of the tetrad e", and is the spin 
connection, which is defined as = |[7 a , ~t b ]e v a ebv^- For the Schwarzschild black hole 

ds 2 = -fdt 2 + jdr 2 + r 2 {d0 2 + siifOdtf 2 ), (2.2) 



with 



r 



where the parameters M represents the mass of the black hole, we can take the tetrad as 



z«=diag{J],±j,r,rsm6), (2.4) 



Introducing an ansatz 



<t>U0M\ , 25) 



with 



/ j+m ym-l/2 v 



j+l-ro ym-1/2 \ 

2 i+ 2 ' for i=Z-- 

2j+2 r / / 

Cho [24] find that the cases for (+) and (— ) in the functions and can be put together, and then the decoupled 
equations can be written as 

d 2 d 2 \ 

dp-W- Vi ) Fi{K > t) = ° ! 1=1,2 (2 - 6) 



where 



dr* 

with 



f» = r + 2Mln (— + —tan^ 1 ( ^-H (2.7) 

V2M / 2w Vl fc l7 

V 1)2 =±^ + W 2 , (2.8) 



r 2 (fc 2 + ^ 2 r 2 ) + ^fcA/2w' 1 ' 
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where A = r(r — 2M). Here k goes over all positive and negative integers. Positive integers represent the (+) case 
of Eq. I|2.8[) with k = j + 1/2 and j = I + 1/2, while negative integers represent the (— ) case of Eq. 1)2. 8[) with 
k = —(j + 1/2) and j = I — 1/2. Equation (|2.8(l shows that the potentials V\ and V2 are related to the metric 
function y 7 / which differs from potentials for the scalar, electromagnetic and gravitational fields [2M |26|. We will 
use the potential l|2.8[) to study the late-time tail behavior of the massive Dirac fields in the Schwarzschild black-hole 
geometry. 

It is easy to see that near the event horizon the effective potential V\ t i becomes zero. Then the radial asymptotic 
solutions can be expressed as 



R, 



(2.10) 



where r* = J 1/fdr is the usual tortoise coordinate. By comparing Eq. 12.1l)fl with the scalar field solutions e ±lujr * 

Q we know that the Dirac field possesses a additional factor e ±5 * a ™ (w). The appearance of the factor seems due 
to we used the "tortoise" coordinate l|2.7|l . However, the factor still exists even if we use the usual tortoise coordinate 
(see appendix for detail). 



III. THE BLACK HOLE GREEN'S FUNCTION 

It is well known that the time evolution of a wave field ^(r,, t) follows from 

*(r„, t) = J [G(r*, K; t)dt*K, 0) + ftG(r., 0)]<K, 
where the black-hole (retarded) Green's function G(r»,r^;i) is defined by 

G(r.y.;t) = 8(t)6{r.-ri). 
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(3.1) 



(3.2) 



The causality condition gives us the initial condition G(r«, r^; t) — for t < 0. In order to get G(r*, rj,; f) we use the 
Fourier transform 



(3.3) 



The Fourier transform is analytic in the upper half lu— plane, and the corresponding inversion formula is given by 



G(n,ri;t) = — 



00 -He 



G(r»,<;c i ;)e^ t d W , 



(3.4) 
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where c is some positive constant. 

We define two auxiliary functions ^i(r*,o;) and ^^(f*,^) which are (linearly independent) solutions to the homo- 
geneous equation 



+ u 2 - V{r) 



%(n,co) = 0, i=l,2. 



Let the Wronskian be 

and using the two solutions and ^2, the black-hole Green's function can be constructed as 



W(uj) \ *i«,w)* 2 (r,,w), r* > <. 



(3.5) 



(3.6) 



(3.7) 
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IV. LATE-TIME TAILS OF THE MASSIVE DIRAC FIELDS 



It is well known that massive tails exist even in a flat spacetime [27|. This phenomenon is related to the fact that 
different frequencies forming a massive wave packet have different phase velocities. One will see that at intermediate 
times the backscattering from asymptotically far regions (which dominates the tails of massless fields) is negligible 
compared to the flat spacetime massive tails that appear here. Hod, Piran, and Leaver [2^ argued that the 
asymptotic massive tail is associated with the existence of a branch cut (in placed along the interval — Li < ui < fi. 
This tail arises from the integral of the Green function G(r*, r'^ui) around the branch (denoted by G c (r*, r'^oj)). So 
our goal is to evaluate G c (r*, r^; u>). 

Now we assume that both the observer and the initial data are situated far away from the black hole so that r ^S> M. 
One may expand the wave-equation (|3.5|) with the potential V\ which is given by Ij2.8|l (or V2) for the massive Dirac 

field as a power series in M/r (neglecting terms of order O ^(M/r) 2 ^) as follows 



with 



dr 2 



£ = (1 - 2M/r) 



4Mw 2 - 2Mll 2 k 2 - ku/fi 



£(r,w) = 0, 



1/2 



1 



1-^ 



»\k\ 



r J 2uo(k 2 + Li 2 r 2 ) 



-1/2 



(4.1) 



(4.2) 



At great distances, the only difference between effective potentials of the scalar and Dirac fields is that the term 
1(1 + l)/r 2 for the massive scalar field is now replaced by (fc 2 — kujj \x)jr 2 for the massive Dirac field. As we will see 
later, it is the difference that leads to new result. 
Let us introduce 



b 2 



2\J \x 2 — uj 2 r — 2wr 

e - Z /2 z l/2+fc $j 

1 , 9 ku 
4 fi 
Mix 2 



- 2Mm. 



Then Eq. I|4.1|l becomes the confluent hypergeometric equation 

d 2 $ 



dz 2 



(l + 26-*)^-(- + 6-o)$ = 0. 



It follows that two basic solutions required in order to build the Green's function can be expressed as 

$1 = Ae" mr (2zur) 1/2+b M(l/2 + 6 - a, 1 + 26, 2zur), 
i> 2 = Be-™ r (2vjr) 1/2+b U{l/2 + b - a, 1 + 26, 2wr), 



(4.3) 
(4.4) 

(4.5) 
(4.6) 



(4.7) 

(4.8) 
(4.9) 



where A and B are normalization constants. The functions M(a,b,z) and U(a, b, z) represent the two standard 
solutions to the confluent hypergeometric equation [28j . U(a, 6, z) is a many- valued function, i.e., there is a cut in ^2- 
Using Eq. I|3.4fl . one finds that the branch cut contribution to the Green's function is described by 



G C (r„r:,t) 



Using the following relations 

^i(2n7re l7r ) 
tp2(2mr) 



1 

2^ 
1 

2~~ 



-p 



IU(we lTr ) 



W{uj) 



— / F{vj)e- lu,t dw. 



(4.10) 



= £ 



> 4 e (l/2+^ e -wr( 2n7T .)V2+«'Af (1/2 + 6 + a, 1 + 26, 2vor) 

r(-26) _ 



r(l/2-6-a) 

£(26) 
r(l/2 + 6-a)' 



(2OT) 1/2+h A/(l/2 + 6 - a, 1 + 26, 2ror) 
(2tur) 1/2 - b M(l/2 - 6 - a, 1 - 26, 2wr), 



(4.11) 
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we find 



W{we i 



-W(w) = AB 



T(26) 
r(l/2 + 6-a) 



Abw. 



(4.12) 



and 



F{w) 



26 



M (i + 6 + a, 1 + 26, 2ot-'JM(~ - 6 + a, 1 - 26, 2<ot%) 



1 



6 - a, 1 + 26, 2ot-*)M( 6 - a, 1 - 26, 2ror 



r(-26)r(^+6-a) 
r(26)r(i-6-a) 



(4tJ7 2 r*<)t+ 6 e 



+b -ro(r*+r^) 



4tu6 



M (- + 6 - a, 1 + 26, 2wr'jM(- + b - a, 1 + 26, 2mr*) 



3 (i+2&)«r M ( + 6 + a, 1 + 26, 2ro<)M(- + b + a, 1 + 26, 2n7r*) 



(4.13) 



If we further assume that both the observer and the initial data are situated in the region M <C r -C M/ (pM) 2 and 
we are interested in the intermediate asymptotic behavior of the Dirac field [M < r < t < M/(M/i) 2 ], we find that 
the frequency range w — 0(y / Ji/t), which gives the dominant contribution to the integral, implies a -C 1. Eq. 14. 1|) 
shows that a originates from the 1/r term which describes the effect of backscattering off the spacetime curvature. 
Thus, the backscattering off the curvature from the asymptotically regions is negligible for the case a <C 1. So, we get 



F{w) 



T{-2b)T{\ + 6) (2vjr,) 1 / 2 + b {2wrj) 1 / 2 + b e-^ r ' + <'> 

' r(26)r(i-6) 4^ 



(1 



,{i+2b)i*j M ( 1 + b,l + 26, 2wr' m )M(- + 6, 1 + 26, 2ror* 



(4.14) 



Noticing that M (a, 6, z) w 1 as z 



F{w) 



where we used the formulae T(2z) 
into the Eq. (|4.10|l . we obtain 



0, we have 

i + e (i+2b)^ r(-26)r(^ + 6) 
r(26)r(i - 6) 



4w6 

7T 1 



(2 W ) 1+2b (r*0 1/2+b 



r(r*r 



sin(7rb) 2 1 + 2b b 2 T(6) : 
i=2 2z - 1 / 2 r(z)r(l/2 + z) and T{ 



I U/2+6 



sin(7T2:) zr(2) 



(4.15) 

Substituting Eq. (|4~T5|) 



G C (r„<,t) = - 



1 1 



= (l + 2fc)i7T 



(r*r 



' U/2+6 



_„ sm(7r6) 2 2b 6 2 



r(6) 2 



w 2 ) h e-^*dw. 



(4.16) 



Unfortunately, the integral can not be evaluated analytically since the parameter 6 = \/l/4 + k 2 — kui/ fi depends 
on to. However, we can work out the integral numerically and the results are presented in figures Hl3l with the form 
In |G c (r*,r», t)\ versus t. Figure ^ describes different k with the same mass (/i = 0.01), which shows that the dumping 
exponent depends on the multiple number of the wave mode. Figure |21 gives different mass /i with the same multiple 
number (k — 1), which indicates that the dumping exponent depends on the mass \i of the Dirac fields. In Fig. |3| 
the lines (a) represent the result of the Green function (|4.16|l of the massive Dirac field with different k, and lines (b) 
show the corresponding result of the massive scalar field. We learn from Fig. 13 that late-time behavior of massive 
Dirac fields is dominated by a decaying tail without any oscillation, and the decay of the massive Dirac field is slower 
than that of the massive scalar field. 



V. SUMMARY AND DISCUSSIONS 



We have studied the intermediate late-time tail behavior of massive Dirac fields in the Schwarzschild black-hole 
geometry and the result is compared with that of the massive scalar fields. We know from the figures that there are 
three differences between the massive Dirac and scalar fields due to the fact that the parameter 6 = y/\/A + k 2 — kuj/fi 
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is related to oj for the massive Dirac fields while b = I + 1/2 for the massive scalar fields, (I) the late-time behavior 
of massive Dirac fields is dominated by an inverse power-law decaying tail without any oscillation, but the massive 
scalar field by an oscillatory decaying tail, (11) the dumping exponent (i.e., a in term l/t a ) for the massive Dirac field 
depends not only on the multiple number of the wave mode but also on the mass of the Dirac field, while a = (I + 3/2) 
for the massive scalar field depends on the multiple number only, and (III) the decay of the massive Dirac field is 
slower than that of the massive scalar field. 
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APPENDIX A: MAY THE CONCLUSIONS CHANGE IN THE USUAL TORTOISE COORDINATES? 

The "tortoise" coordinate (|2.7(l is a function of the background geometry and the test field, which is different from 
the usual tortoise one, r* = J Might the conclusions of this paper change if we use usual tortoise coordinate 
instead of the coordinate (|2.7|l ? In the appendix we will address this question. 

The Dirac equations are given by |29| 

V2V BB ,P B +ipQ B , -0, 

y/2V BB ,Q B + iiJLP B ,=0, (Al) 
where Vbs' is covariant differentiation, A BB / is the electromagnetic field potential, P B and Q B are the two-component 



pot i 

spinors, and [i is the particle mass. In the Newman-Penrose formalism |30j the equations become 

(D + e- p)P Q + (5 + 7T - a)P x = 2- 1/2 ifiQ 1 ', 
(A + fi- 7 )P X + {5 + /3- t)P° = -2-^0°', 
(D + e-p)Q ' + (6 + n- a)Q 1# = -2- 1/2 ipP\ 

(A + fi - 7)Q X ' + (5 +fi - f)Q ' = 2- 1 / VP°, (A2) 
For the Schwarzschild spacetime the null tetrad can be taken as 

i" = (^, o, 0, 0), 



1 A 

2 

1 



(1, 0, 0) 



y/2r V ' ' s i n ® 
Then, if we set 

r 

Q 1 ' = g 1 {r,6)e-' l{uJt - mip) 1 



0, 0, 1, — , (A3) 



0°' = --g 2 [r,e)e- i ^ t - mv \ (A4) 



r 



7 



where u> and m are the energy and angular momentum of the Dirac particles, Eq. i|A2|) can be simplified as 



£>o/i 



l 



V2 



A/2/2 



1 



v/2 



(A5) 



with 



d iK r - M 
d iK r-M 
d rr 



89 sin( 
d m 



n cot I 



n cot I 



50 sin 6» 
X = r 2 oj. 

It is now apparent that the variables can be separated by the substitutions 

fi = R- 1/2 {r)S_ 1/2 (9), 
f2 = R + i /2 (r)S +1/2 (0), 

91 = R+i/2(r)S^ 1/2 {0), 

92 = R-i/2{r)S +1/2 {9). 



(A6) 



Thus, we have pi 



-R-1/2 — 0, 



(A7) 



(A8) 



1 / 2 ~" u A + i/xr" 

and \fKR + i/ 2 satisfies the complex-conjugate equation. The decoupled equations can then be explicitly expressed as 

-dR_ 1 / 2 s 



dr 
1 d 
A dr 



A 3/2 



dr 

dR+1/2 
dr 



with 



P+ = 



K 2 + i(r — M)K 
A 

K 2 - i(r - M)K 



liior 



i/j,A dR_i/2 
A + ifir dr 
ifiA dR +1 / 2 
X — i[ir dr 



P-R- 1/2 = 0, 
+ P+R +1 ,2 = 0, 



A 2 - * 



A + ijir 



^r + 2s + ^ r r M) - K) -, 2 r 2 -A 2 , 



A A — i[ir 

where A 2 = (I — s)(l + s + 1) is the separation constant. Introducing an usual tortoise coordinate 



r A dr > 



(A9) 
(A10) 

(All) 
(A12) 

(A13) 



and resolving Eqs. I|A9|I and l|A10|) in the form 



R+i 



A-V4 



/2 



R-i 



r 

A 1 / 4 



(X 2 + f i 2 r 2 ) 1 / 4 e- M ^, 



/2 



(X 2 +LL 2 r 2 ) 1 ' i e M fH-. / 



(A14) 
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with 



?9 = arctan(/xr/A), 



we obtain two wave-equations 



drl 



+ 



d 2 *_ h f di7 - 
dr 2 \ dr* 



where 



H- 



H% + -P+ }V+ = 0, 



H 2 _ + !> = 0, 



if i 



1 dA A 

4r 2 dr r 3 2(A + i/ir) r 2 ' 
1 dA A i/z A 



4r 2 dr r 3 2(A — ifir) 



Near the event horizon the asymptotic solutions are 



(A15) 

(A16) 
(A17) 

(A18) 
(A19) 

(A20) 



One may expand the wave-equations i|A16(l and i|A17l) as a power series in M/r (neglecting terms of order O ^(M/r) 2 ^ ) 
as follows 



2 i A, 



dr 2 



Ijj - (J, + 



AMuj 2 - 2Mfi 2 A 2 + 



£fc =0, 



(A21) 



where ^± = (1 — 2il//r) 1 ' 2 f ±. It is of interesting to note that Eq. I|A21(1 becomes Eq. I|4.1|l if we replace A by 
— k, which shows that the conclusions of this paper do not change if we use usual tortoise coordinate instead of the 
coordinate (j2.7(l . 
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FIG. 1: The figure describes In |G c (r*, r'„, t)\ versus t for different k with the same mass (fj, = 0.01), which shows that the 
dumping exponent depends on the multiple number of the wave mode. 




FIG. 2: The figure describes In |G c (r*, r'„ t)\ versus t for different mass /i with the same multiple number (k = 1), which 
indicates that the dumping exponent depends on the mass fi of the Dirac fields. 



11 




500 1000 1500 2000 2500 t 500 1000 1500 2000 2500 t 




500 1000 1500 2000 2500 t 500 1000 1500 2000 2500 t 




1000 2000 3000 4000 t 1000 2000 3000 4000 t 



FIG. 3: The above four figures describe In |G c (r*, r'„, t)\ versus t for fi = 0.02, and below four for fi = 0.01. The dashing 
lines represent the result of the Green function of the massive Dirac field with different k. For comparing, we also show the 
corresponding result of the Scalar field with solid lines. It is shown that the late-time behavior of massive Dirac fields is 
dominated by a decaying tail without any oscillation, and the decay of the massive Dirac field is slower than that of the massive 
scalar field. 



